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ABSTRACT 

(N 

The announcement of 2300+ candidate transiting exoplanets (KOIs) orbiting ~1800 

host stars discovered with the Kepler mission enables a plethora of ensemble analysis of 



3 



the architecture and properties of exoplanetary systems. We use the transit durations 
to probe the ensemble validity of stellar parameters for Kepler candidate host stars. 

q^ While the new KOI stellar parameters are an improvement over the previous release, 

our analysis shows that a systematic trend in the ensemble stellar radii remains, which 

n/ affects the KOI exoplanet radii inferred from transit depths. We also compare the distri- 

bution of observed transit durations to a model generated from the known eccentricities 
of radial velocity (RV) discovered exoplanets. We find that the Kepler and RV distri- 
butions differ at a statistically significant level, even after accounting for the transit 

probability as a function of the eccentricity, inclination, and periastron angle. Addi- 
.£3 

tionally, smaller radii KOIs are preferentially found on more circular orbits. Finally, 

there is an over-abundance of KOIs with long transit durations. We identify that errors 

in stellar metallicity and possibly extinction, resulting in under-estimated stellar radii, 

£> can account for these sources rather than eccentric planets improbably transiting near 

oo 
oo 

en 

o 



a 



apastron. 

Subject headings: exoplanets - statistical 



1. Introduction 



The Kepler mission is revolutionizing our understanding of exoplanets (Borucki et al. 2010), 



including among its many existing and assuredly forthcoming highlights the discovery of three 



terrestrial exoplanets orbiting the M dwarf Kepler Object of Interest (KOI) 961 (Muirhead et al 



2012), and a 2.4 R§ exoplanet in the habitable zone of its host star (Borucki et al. 2012). A list 



of 312 KOIs was published in Borucki et al. (2011a), derived from Q0-1 Kepler time-series data. 
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This was soon followed by a second list of 1235 KOIs in Borucki et al. (2011b) derived from Ql 



Q5 Kepler data. The first two KOI releases relied on stellar parameters from the Kepler Input 



Catalog (KIC, |Brown et al.||201l| ). |Batalha et aT] ( |2012[ ) (hereafter B12) recently announced 2321 
candidate transiting exoplanet KOIs orbiting 1783 host stars from an improved pipeline analysis 
of the Q1-Q6 Kepler data, and provided updated stellar host parameters. This large ensemble of 
exoplanet candidates, including many multi-exoplanet systems, beckons for the ensemble analysis 
of exoplanet ary system architectures ( Fabrycky et al.||2012 Howard et al.||2011 Morton & Johnson 



2011; Plavchan k Bilinski 2012; Figueira et al. 2012 Youdin 2011). 



The ensemble analysis of KOIs and the determination of the frequency of Earth-sized planets 
(77®) in habitable zone orbits relies on the accuracy of the estimated physical stellar host parameters 
such as mass, radius and effective temperature (Batalha et al. 2010a|b Pinsonneault et al. 2012 



Traub 2012). Recent studies suggest that the Kepler Input Catalog (KIC) mischaracterizes some 
objects, in particular the stellar surface gravity for stars with effective temperatures less than 



4500 K. For example, Mann et al. (2012); Ciardi et al. (2011) demonstrate that M "dwarf" KIC 



targets brighter than a Kepler magnitude of 14th are most likely giants, and Muirhead et al. (2011) 



reported the over-estimate of M dwarf radii for KOI host stars from the second KOI candidate 



release in Borucki et al. (2011b). B12 acknowledges these limitations of the new and old stellar 
parameters. 

In this work we use the transit durations of KOIs for two purposes: one, as a probe of the 
inferred eccentricity distribution of the candidate exoplanets in comparison to exoplanets discovered 
with the radial velocity (RV) technique, and two, as a new diagnostic of the accuracy of the 
estimated stellar host parameters. The analytic framework for these two investigations is already 



laid out in Ford et al. (2008) (hereafter F08) and Burke (2008) and references therein. A more 



comprehensive analysis is also already carried out with the Borucki et al. (2011b) KOI list by 



Moorhead et al. (2011). Additionally, Dawson & Johnson (2012) carry out a thorough analysis 



when the impact parameter can reliably be estimated for Jovian-sized exoplanet KOIs with high 
S/N transit time-series. In particular, we confirm that the new KOI stellar parameters and KIC 
values still have some systematic errors. We also find that the distribution of the KOI transit 
durations is in disagreement with transit durations modeled from the eccentricity distribution of 
RV-discovered exoplanets. 



2. Samples 

We use the second and third tabulation of KOI candidate exoplanet and stellar host parameters 



from Borucki et al. (2011b) and B12 respectively, the latter of which is currently available at the 



NASA Exoplanet Archive (http://exoplanetarchive.caltech.edu/). Details on the methods used to 
tabulate these parameters can be found in those references and citations therein. We assume that 
there are no false-positive KOIs, and for simplicity we treat multiple KOIs orbiting the same host 
star as independent statistical tests. We make additional subsets of terrestrial, Neptune-like, and 



Jovian KOIs at exoplanet radii of R p i <2, 2< R p i <6, and R p i >6 Rq earth radii respectively. 
We do not filter the KOIs by orbital period or by their host mass, temperature, or radius. We do 
exclude KOIs from our analysis where only one transit is observed and the tabulated period is only 
a lower bound. 

We also use the eccentricity values of 428 published and confirmed RV-discovered exoplanets 
as listed at the NASA Exoplanet Archive as of March 7th, 2012, with a mean eccentricity of 0.25. 
For simplicity, we do not filter the RV exoplanets by planet mass, period, or stellar properties, 
but we note that both samples are reasonably complete for large planets up to periods of ~8 
months. Whereas the distribution of periastron angles for RV exoplanets is uniform random on 
a sphere, the distribution of periastron angles for transiting exoplanets will more strongly favor 
angles perpendicular to the plane of the sky for eccentric planets due to the increased probability 



of transit (F08, Burke 2008). Thus, we do not use the measured periastron angles for RV planets 



since they have an incorrect distribution for comparison to the transiting planets. 



3. Observed Transit Duration Anomalies 

Equal to tq in F08 in the limit R p i/R* — > 0, we introduce the transit duration anomaly 
dimensionless parameter a: 

_ aTir 

° ~ (R* + R pl )P (1) 

where a is the ratio of the observed transit duration to the expected transit duration for a circular 
edge-on orbit, and can be computed from the KOI table parameters in B12. R* and R p i are the 
stellar and exoplanet radii respectively, P is the orbital period of the exoplanet, T is the duration 
of the transit, and a is the semi-major axis of the exoplanet. The semi-major axis of the exoplanet 
a is in turn calculated from the stellar mass M*, the orbital period of the exoplanet P, and Kepler's 
third law. P and T are observed quantities directly measured from the time-series, whereas R* and 
M* (and consequently a) are estimated from ancillary observations and models such as observed 
colors, magnitudes, spectroscopy and fits to theoretical isochrones (Batalha et al. 2010a|b 2012 



Muirhead et al. 2011). Values for the transit duration anomaly a are tabulated in Table 1. 



For an eccentric, non-edge-on orbit, the square of the transit duration anomaly can be written 
in a similar fashion to the square of Equation 1 in F08 as: 



2 _ I d t\l I 
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where e is the eccentricity, R p is the exoplanet radius, b is the impact parameter, and dt is the 
exoplanet-star separation during transit as defined in F08. Assuming as in F08 that a/R* ^> 1, 



a/Rp 3> 1, and T/P <C 1, dt can be assumed to be approximately constant during transit (hence 
the approximation rather than equality in Equation 2). 

From Equation 2, the eccentricity distribution of transiting exoplanets can be inferred from 
the distribution of observed transit duration anomalies (e.g. the mean, standard deviation), as 
discussed in great detail in previous works in Burke ( 2008| ); Ford et al. (2008); Moorhead et al 



(2011); Barnes (2007). a is the correct quantity to use as a proxy for eccentricity instead of the 
transit duration excess in hours, since transit durations increase with semi-major axis and period 
due to slower orbital velocities at a fixed eccentricity. 

In particular, we can identify potential eccentric exoplanets with values of a <C 1 transiting 
near periastron ora>l transiting near apastron. Since b is degenerate with eccentricity, the KOIs 
with a <C 1 may alternatively be grazing transits. We tabulate these candidates in Table 2. We 
find in §5 that most KOIs with q>1 possess errors in their stellar parameters, and are excluded 
from Table 2 unless log[Fe]/[H] > -0.11 and A v < 0.33 mag. 

The four KOIs with the smallest inferred transit duration anomalies compared to a circular 
orbit - KOIs 338.01, 338.02, 977.01, and 1054.01 - are reported to orbit giant stars with log g < 2.5 
with periods of 1.4-7 days interior to the estimated stellar radii, and may instead be associated with 
photospheric activity or a blend. Visual inspection of the time-series for KOIs 977 and 1054 support 
this hypothesis. The KOI 338 system has two candidate exoplanets, exhibits clear transit-like dips 



that are short w/r/t to the transit period, and is specifically mentioned in B12 and Fabrycky et al 



(2012). We also exclude these four KOIs from Table 2. 



The mean value for a for all KOIs from B12 is 0.921. The expected value for a assuming all 
exoplanets are on circular orbits with a random impact parameter between and 1 with < b >= 0.5 
is v3/2 ~ 0.866. The higher and statistically significant observed mean value of a suggests that 
there exists a systematic ensemble under- (over-) estimate of stellar radii (mass). 

To ascertain the validity of this claim, we first briefly discuss the relevance of detection biases 
in the list of KOIs. First, one might expect a deficiency of transits at high impact parameters 
(i.e., < b > less than 0.5) due to the lower S/N for shorter transit durations, and due to the 
removal of V-shaped grazing transits from the list of KOIs because of increased probability of 
false-positives from stellar eclipsing binaries. B12 includes estimates of the impact parameter b, 
which can be determined from the ingress and egress times for high S/N transits after accounting 



for limb- darkening (Seager & Mallen-Ornelas 2003). This detection bias is only readily apparent 



for KOIs with b > 0.9 as tabulated in B12, and it is not correlated with exoplanet radius. In fact, 
because a number of the KOIs are eccentric and transiting near periastron (as opposed to circular 
orbits), the distribution of calculated impact parameter in B12 has a median value of 0.7, which is 
greater than 0.5. This larger value for < b > suggests both that Kepler can adequately detect short 
duration transits, and that some of the KOIs must be on eccentric orbits, further exasperating the 
discrepancy between the observed mean value of a for all KOIs and that expected for a population 
of exoplanets on circular orbits with uniform random impact parameters. We do not include the 



effects of limb-darkening on our analysis. We also do not use the impact parameters in B12 since 
a circular orbit is assumed, and there are large uncertainties in constraining b due to the Kepler 



30-minute cadence (Fabrycky 2012). Regardless, we find that adjusting the upper bound to the 



assumed impact parameter distribution in our analysis does not alter our results. 

Additionally, one might expect that there is a deficiency of small exoplanets with short transit 
durations relative to larger exoplanets and/or longer transit durations at higher S/N, leading to 
a bias towards larger a values for smaller exoplanets. However, this bias does not appear to be 
present for the KOI parameters tabulated in B12. All but 7 KOIs have a transit detection S/N>7, 
and there is no observed trend of transit duration in hours as a function of planet radius. This lack 
of a bias in our analysis is due to the long Kepler baseline, and the fact that transits are longer 
for longer orbital periods. However, we do not draw any conclusions regarding the frequency or 
distribution of exoplanets as a function of orbital period, which would necessarily have to take this 
bias into account. 



4. Comparison to Radial Velocity Discovered Exoplanet Eccentricities 

Since the impact parameters and angles of periastron for KOIs are generally not constrained, 
Equation 2 does not permit a direct quantitative statistical comparison to the eccentricity distribu- 



tion of RV discovered exoplanets. As recognized by Burke (2008), F08 and others, the distribution 
of transit duration anomalies a can be used as a proxy instead. We convert the distribution 
of eccentricities of confirmed RV exoplanets into a simulated distribution of a transit duration 
anomalies. We assume that every radial velocity exoplanet "transits" with a distribution of impact 
parameters b and angles of periastron u that are properly weighted by the transit probability. As 
noted previously, the angle of periastron cannot be assumed to be uniform random on a sphere as is 
true for RV discovered exoplanets. An accurate probability distribution for the impact parameter, 



eccentricity, and angle of periastron for transiting planets is given in Burke (2008, Equations 14-16; 
Figure 4). We numerically integrate the transit probability using the prescription in (Burke 2008) 
combined with the known RV sample eccentricity values to derive a simulated distribution of a 
values for the RV sample. We then directly compare this to the CDFs of a values for all KOIs, 
KOIs with R p i < 2i? e , R® < R p[ < R®, and R p i > 6i? e . We plot these CDFs in Figure 1. 

The simulated RV sample a distribution can be approximated by a normal distribution with 
a mean and standard deviation of (/j,, a) = (0.714,0.305) (x 2 =0.659). This result is consistent 
with the models in F08 for an assumed Rayleigh distribution of eccentricities with an average 
eccentricity of ~0.2-0.3. This also verifies that we correctly computed the a distribution from 
the eccentricities. A large value of a > 1 (assuming accurate stellar parameters) requires a high 
eccentricity and transit near apastron, which is a low probability occurrence. This explains the 



notable lack of simulated transit duration anomalies with a >1.5, which is consistent with Burke 



(2008 Figure 4). 



The a distribution for all KOIs with (ft, a) = (0.921,0.288) (x 2 =2.03) is significantly different 



from the RV a distribution, in disagreement with Kane et al. (2012). We calculate distribution mean 
and standard deviations of (p,<r) = (0.978,0.255) for R pl < 2i? e , (0.901,0.279) for 2R (B < R P i < 
6fi©, and (0.761, 0.331) for R pl > 6.R© (x 2 =0. 669, 1.02,0. 556). The differences between almost all of 
these distributions are statistically significant. Formally, we evaluate the 1-D Kolmogorov-Smirnov 
(K-S) test statistic between all distributions. The a distribution for RVs has a probability of being 
drawn from the same parent population as the KOIs with R p i > 6i?e (2Rq < R p i < 6R^,R p i < 
2i? e ,all) of 4.2 x 10~ 4 (4.0 x 10~ 23 ,4.7 x 10- 40 ,1.5 x 10" 30 ). The KOIs with 2i? e < Rpl < ®R(B 
and Rpi < 2R^ have a probability of 13% of being drawn from the same parent population (< 2 a) 
and are thus the most similar, whereas the KOIs with R p i > 6i?® are statistically distinct from 
the KOIs with 2i? e < R p i < 6i? e and R p i < 2i? e with K-S test probabilities of 3.0 x 10~ 9 and 
1.4 x 10 -19 respectively (> 5 a). In other words, KOIs with R p i > 6i?© most closely resemble the 
RV distribution perhaps as expected, but are still statistically distinct at the > 3 a level (and with 
higher statistical significance for the remaining KOIs). This is also visually apparent from Figure 
1. Similarly, the KOIs with R p i < 6R^ appear to be drawn from a distinct eccentricity distribution 
than KOIs with R p i > 6-R®, with little difference, if any, above and below R p i = 2R®. 

From this we can conclude that smaller radius KOIs are increasingly more likely to have 
smaller average eccentricities within these samples. This apparent trend may solely be due to the 
incompleteness of Kepler at long orbital periods for smaller exoplanets (§3) combined with the well- 



established tidal circularization of exoplanets at short orbital periods (Butler et al. 1997 Rasio & 



Ford 1996). Alternatively, Fabrycky et al. (2012) notes that smaller exoplanets are more likely to 



be found in multiple exoplanet co-aligned systems; the smaller relative mutual inclinations implies 
smaller average eccentricities. 

Finally, the a < 0.75 distribution of KOIs with R p i > 6R& most closely (and encouragingly) 
resembles the RV exoplanet distribution, as might be expected. However, there is a higher occur- 
rence rate of KOIs with a >1.5 than expected for both a Rayleigh eccentricity distribution (F08) 
and for our simulated RV sample distribution. This appears to be especially true for Jovian KOIs 



with R p i > 6i?0. This excess was also noted in Moor head et al. (2011) and Wang & Ford (2011) and 
we confirm that this excess is still present with the revised stellar parameters and new list of KOIs. 
If we assume that our RV and KOI samples are drawn from the same population of exoplanets, 
some of these KOIs must be attributable to errors in the host stellar parameters or alternatively 
are false-positives. It is not clear how the stellar radius (mass) could be under- (over-) estimated 
by factors exceeding 1.5 (1-5 3 ) to correct some of these a-values. However, a blanket rejection of 
these KOIs is not recommended since that could inadvertently exclude a genuine highly eccentric 
exoplanet transiting near apastron. 



5. Accuracy of Stellar Parameters Inferred from Transit Duration 

The significant difference between the transit duration anomaly a distributions for RV exoplan- 
ets and Kepler candidates raises the question about the validity of the KOI host stellar parameters. 
In Figure 2 we plot the dimensionless a parameter as a function of stellar mass, radius and effective 
temperature, for both the second and current lists of KOIs. We bin the a values for each stellar 
parameter, with bins of 250 K in effective temperature, 0.1 Mq in stellar mass, and 0.1 Rq in 
stellar radius. Our results are insensitive to adjustments in the bin width. We derive median values 
and quartile ranges for each bin, which are over-plotted in Figure 2 and show distinct trends as a 
function of stellar mass and radius for both the old and new KOI lists. To quantify these trends, 
we fit a line to the median values, excluding bins with fewer than 4 KOIs. 

The new list of KOIs yields less scatter in the median values of a as a function of stellar 
parameters, indicating that the revised stellar parameters in B12 are indeed improved. However, 
we identify a statistically significant trend for decreasing stellar radius (mass) of a 0.27 (0.28) change 
in a per Rq (Mq) at the ~6-er (3-<r) level for both the current and old KOI lists. The ensemble 
a values are not found to be linearly biased as a function of stellar effective temperature. The 



identified trends in mass and radius may have been overlooked in Moorhead et al. (2011) because 
they only examined the CDFs of the a distribution as a function of T e ff. Additionally, the linear 
trend in stellar temperature for the old KOI list is statistically significant if one excludes KOIs with 
T e ff > 7000 K where the median a values decrease again (not shown in Figure 2). This trend is 
not seen for the new KOI stellar effective temperatures. 

The trends in stellar mass and radius imply one of two scenarios. First, there could be a 
systematic error in these two parameters for the ensemble of KOIs. If we assume that the median 
a value for 1.0 Rq (Mq) KOIs is accurate, then the ensemble radii (or masses) are under- (over- 
estimated by ~15% (~50%) at 0.5 Rq (Mq), and the ensemble radii (or masses) are over- (under- 
estimated by ~15% (~50%) at 1.5 Rq (Mq). The larger implied error in stellar mass is due to 
the weaker 1/3-rd power dependence of a on stellar mass. The ensemble stellar masses in B12 and 
the KIC are likely not in error by ~50% at the low- and high- mass ends. The a-trend bias is thus 
largely due to systematic ensemble errors in the stellar radius at the low- and high-radius ends by 
the indicated amounts. 

The first scenario relies on the assumption that the eccentricity distribution of exoplanets is 



intrinsically stellar-mass independent. Alternatively, Howard et al. (2011) suggests that terrestrial 
planets are more common around lower-mass stars, and Figure 1 demonstrates that terrestrial 
planets may preferentially be found on more circular orbits. This picture would at first appear to 
be consistent with the observed trends in a both as a function of stellar mass and radius. In other 
words, the systematic "error" in a could instead be due to a real change in eccentricity distributions. 
A ~15% change in a between 0.5 and 1.0 M (or 1.0 and 1.5 Mq) would correspond to a >0.2 
change in the average ensemble eccentricity of e ~(0,0.25, 0.6) for (0.5,1,1.5) Mq as inferred from 
Figure 9 in F08. Such a large trend in average eccentricity is not observed for the RV exoplanets 



in our sample. 

Unfortunately, these two scenarios are degenerate, so we cannot rule out a modest change in 
ensemble eccentricity distributions as a function of stellar mass. However, given the lack of a large 
eccentricity dependence on stellar mass for RV discovered exoplanets, we conclude that the bulk 
trends in ensemble a as a function of stellar mass and radius for KOIs are likely entirely due to 
systematic errors in stellar radii. The ensemble KOI under-estimate of stellar radii for K/M dwarf 



runs contradictory to the results of Muirhead et al. (2012). 



Finally, the over- abundance of a values greater than 1.5 identified in §4 is independent of stellar 
mass, radius and temperature. We compare KOI a values against every property calculated in the 
KOI tables in B12. We find that KOIs with a > 1.5 are preferentially found around low metallicity 
KOIs with log[-Fe]/[iI] <-0.11 for all planet radii, and for Ay > 0.33 mag for R p i > 6i?© as shown 
in Figure 3. The cutoff values for log [Fe] /[-£/] and Ay correspond to the sample medians. While 
the extinction may be a crude proxy for brightness, there is no a priori reason to expect these two 
trends. Thus we conclude that the over- abundance of KOIs with a > 1.5 is due to errors in the 
calculated stellar metallicities (and possibly extinction) for these sources, rather than constituting a 
genuine class of exoplanets with high eccentricities (improbably) transiting far outside of periastron. 
These errors in turn imply under-estimates of the stellar radii for these KOIs of >50%. Instead, 
these stellar hosts are likely more distant (sub-)giants with larger stellar radii and larger secondary 
companions. 



6. Conclusions 

We have carried out an updated analysis of the transit duration anomalies with the recently 
released list of new Kepler exoplanet candidates. In particular, we looked at the KOI distribution 
of transit durations compared to what would be expected from the eccentricity distribution of RV 
discovered exoplanets as a function of stellar host parameters. We find systematic errors in the 
stellar radius and metallicity preclude a quantitative ensemble comparison of the two samples at this 
time. However, we find that terrestrial planets may preferentially be found on more circular orbits 
compared to Jovian planets. The systematic bias in stellar radius impacts the inferred distributions 
of exoplanet properties, including radius and habitability. 
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The authors thank John Johnson, Dan Fabrycky, Philip Muirhead, Andrew Youdin, and Alex Hayes 
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Table 1. KOI Transit Duration Anomaly and Eccentricity Limits 



KOI 



1.01 


0.589 


0.485 


2.01 


0.666 


0.386 


3.01 


0.903 


0.102 


4.01 


0.444 


0.671 


5.01 


0.467 


0.642 


5.02 


0.776 


0.248 


7.01 


1.22 




0.01 


0.710 


0.330 


2.01 


1.12 




3.01 


0.658 


0.396 



0.198 
0.117 
Full tabic available online in electronic format. 

a For a transit assumed to occur at pcriastron 
with a zero impact parameter. 

b For a transit assumed to occur at apastron with 
a zero impact parameter. 



Table 2. Candidate Eccentric KOIs 



KOI 


a 


e p a 


p b 


1845.02 


0.155 


0.953 






371.01 


0.156 


0.952 






2519.01 


0.178 


0.939 






2522.01 


0.190 


0.931 






2287.01 


0.199 


0.924 






403.01 


0.206 


0.918 






1815.01 


0.221 


0.907 






1164.01 


3.90 




0.877 


2046.01 


2.47 




0.719 


783.01 


2.41 




0.7077 



a For a transit assumed to occur at 
periastron with a zero impact param- 
eter. 

b For a transit assumed to occur at 
apastron with a zero impact parame- 
ter. 
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Fig. 1. — Cumulative Distribution Functions for: the distribution of transit duration anomalies 
a modeled from the distribution of RV-discovered exoplanet eccentricities (black) with a mean 
and standard deviation of (//, a) = (0.714,0.305), the distribution of a for all KOIs (red) with 
(p,a) = (0.921,0.288), R pl < 2i? e KOIs (blue) with (p,a) = (0.978,0.255), 2,R ffi < R pl < 6^® 
KOIs (green) with (fi,a) = (0.901,0.279), and R pl > 6i? e KOIs (pink) with (fi,a) = (0.761,0.331). 
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Fig. 2. — Distributions for the transit duration anomalies a as a function from left to right of stellar radius, mass, and effective 
temperature for the new KOI list in the top row, and old KOI list in the bottom row, respectively. Overplotted are binned 
median, lower and upper quartile ranges, with bin sizes of 0.1 Rq, 0.1 M©, and 250 K. Some bins are missing upper/lower 
quartile ranges since only a couple stars fall into that particular bin. Linear fits to the binned median values (not shown) yield 
slope coefficients clockwise from top left of -0.27±0.05, -0.28±0.09, -(3.2±2.6)xl0- 5 , -(8.3±4.7)xl0~ 5 , -0.28±0.09, & -0.33±0.05 
respectively with the appropriate units, indicating statistically significant ensemble trends for stellar radius and mass. 
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Fig. 3. — Distributions for the transit duration anomalies a as a function of stellar metallicity on 
the left and extinction on the right. KOIs with exoplanet radii of R p i < 2Rq, 2R^ < R p i < 6Rq, 
and R p i > 6R& are shown in blue, green and red respectively. The horizontal line at a = 1.5 
corresponds to a CDF of ~100% for simulated a values from the known eccentricity distribution 
of RV-discovered exoplanets. This indicates that most a > 1.5 KOIs likely have large systematic 
errors in stellar parameters. The vertical lines correspond to the median values for all KOIs in 
metallicity (-0.11) and extinction (0.33). 



